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investigation using mitochondrial haplotypes 

Manja U. Bohme, Norbert Schneeweiss, Uwe Fritz, Ml Moravec, Igor Majlath, 
Viktoria Majlathova & Thomas U. Berendonk 

Abstract. The range of Lacerta viridis viridis is the result of a rapid postglacial dispersal, followed by 
range-restriction events leading to the isolation of the present northern relict populations. Current 
conservation efforts focus on these endangered, northern isolates. An understanding of the genetic diversity 
of these populations is a prerequisite for any conservation measures. In order to provide such data we 
analysed mtDNA sequence variation in 57 L. v. viridis from the northernmost part of the subspecies range, 
including representatives of isolated relict populations in Brandenburg (Germany) and northern Bohemia 
(Czech Republic). We detected a positive correlation between genetic distances of mtDNA sequences and 
geographic distances, which corresponds well with the migration behaviour of L. v. viridis. Furthermore, 
our data reveal the importance of a comparative analysis for estimating the genetic diversity of an 
endangered taxon like Lacerta v. viridis. 

Key words: Reptilia: Lacerta viridis viridis ; edge population; mitochondrial haplotypes; cytochrome b; 
control region; genetic distance. 


Introduction 

In modem conservation strategies peripheral 
populations become more and more impor¬ 
tant. Although central populations of a spe¬ 
cies are usually not endangered, the more 
threatened peripheral populations represent 
important sources for adaptation and for the 
evolutionary future of a species (Lesica & 
Allendorf 1995). Peripheral populations of¬ 
ten inhabit atypical, less suitable habitats 
compared to populations in the central part 
of the species’ range. Consequently, in pe¬ 
ripheral populations individuals are often 
subject to strong environmental pressure, 
different from the situation in the centre of 
the species’ range. These conditions near the 
range border are thought to trigger the evo¬ 
lution of genetic and/or morphological dis¬ 
tinctiveness (Lesica & Allendorf 1995). 
Their localisation at the range border and 
their smaller population size makes these 
populations vulnerable to founder effects, 


inbreeding and genetic drift. It is well known 
that decreasing genetic variation can lead to 
a higher extinction risk, especially in small 
populations. These problems are the focus of 
many recent conservation studies (Prior et al. 
1997, Mockford et al. 1999, Edenhamn et al. 
2000, Garner et al. 2004). 

The Green Lizard Lacerta viridis (Lau¬ 
renti, 1768), is one species out of eight in the 
genus Lacerta sensu stricto (Nettmann 2001). 
Currently five subspecies of L. viridis are 
recognized (Rykena et al. 2001): L. v. viridis 
(Laurenti, 1768), L. v. guentherpetersi (Ryke¬ 
na et al., 2001), L. v. infrapunctata (Schmidt- 
ler, 1986), L. v. meridionalis (Cyren, 1924) 
and L. v. paphlagonica (Schmidtler, 1986). 
Lacerta v. viridis inhabits a wide range, ex¬ 
tending from the Balkan Peninsula north¬ 
wards across the Carpathian Basin to the 
more isolated edge populations in the Czech 
Republic and the northernmost populations 
in eastern Germany. This subspecies is high¬ 
ly endangered in Germany and the Czech 


© 2006 Deutsche Gesellschaft fiir Herpetologie und Terrarienkunde e.V. (DGHT) 
http://www.salamandra-journal.com 


29 










Manja U. Bohme et al. 


Republic (Elbing 2001a, Mikatova 2001) and 
many conservation efforts are undertaken to 
protect populations. 

The German populations are situated in 
two narrow areas: one near Passau (Bavaria) 
and the other in the Lusatia region (Nieder- 
lausitz, Brandenburg). The Brandenburg ha¬ 
bitats differ significantly from typical L. v. 
viridis habitats. The isolated Lusatian popu¬ 
lations represent relict populations resulting 
from a rapid range expansion from a southern 
glacial refuge during the Holocene warming. 
According to Peters (1970), the migration of 
this mesophile lizard to eastern Germany 
started 8000-7000 years BP in Preboreal 
times. Natural changes such as postglacial 
climatic oscillations, vegetation succession 
and anthropogenic impact on habitats have 
probably caused an increased habitat frag¬ 
mentation. This process is leading to the 
current situation, with patchy localities har¬ 
bouring only small isolated populations as 
in eastern Germany (Brandenburg). Due to its 
rare occurrence in Brandenburg, L. v. viridis 
has provoked attention since the 1930s 
(Hecht 1930, Mertens & Schnurre 1946, 
1949, Elbing 1996, 2001a, b, Schneeweiss et 
al. 2004). These investigations revealed the 
collapse of most of the Brandenburg popula¬ 
tions during the 20th century. A current con¬ 
servation programme in Brandenburg moni¬ 
tors the remaining populations and includes 
a captive breeding group for reintroduction 
measures in regions where the subspecies has 
already become extinct (Schneeweiss 2001). 
Later we will refer to this breeding group as 
a captive population. Because of the endan¬ 
gered situation in Brandenburg, information 
on the phylogeographic origin and the ge¬ 
netic variation of the wild populations be¬ 
comes a necessity. Our investigation on the 
genetics of L. v. viridis is embedded in this 
program and seeks to compare the genetic 
variation of Brandenburg Green Lizards with 
other northern and central populations. 

As a genetic marker we use the fast evolv¬ 
ing mitochondrial DNA (mtDNA) fragment 
containing the cytochrome b gene (cyt b) 
and the control region (CR). 


These DNA fragments have been used as 
powerful tools for the detection of phylogeo¬ 
graphic patterns and Holocene range fluctua¬ 
tions so far (Lenk et al. 1999, Haring et al. 
2000, Seddon et al. 2001, Babik et al. 2004). 
Furthermore, these mitochondrial fragments 
yield information about the genetic relation¬ 
ships of populations and their genetic varia¬ 
tion even within small regions (Gubitz et al. 
2000, Paulo et al. 2002, Hirota et al. 2004). 

To test the genetic status of the endan¬ 
gered L. v. viridis populations in Branden¬ 
burg, we estimated the genetic differentia¬ 
tion of these populations from the geogra¬ 
phically closest populations of the 
continuous range (edge populations) in the 
Czech Republic. We also estimated genetic 
differentiation in comparison to the more 
central populations of the subspecies’ range 
in Slovakia, Austria and Hungary. To obtain 
more information about the genetic status of 
the endangered populations within the Bran¬ 
denburg region, we conducted an analysis of 
the population differentiation between cap¬ 
tive and wild Brandenburg populations. This 
seems to us a prerequisite for further projects 
concerning the reintroduction of L. v. viridis 
in Brandenburg. 


Materials and methods 

Taxon sampling and DNA extraction 

Blood samples from 57 individuals of L. v. 
viridis from the northern and more central 
parts of the species’ range were collected 
(Tab. 1). We sampled three isolated relict 
populations and one captive breeding popu¬ 
lation from eastern Germany. Furthermore, 
four populations from the Czech Republic, 
one population from Slovakia, three popula¬ 
tions from the eastern part of Austria and one 
population from northern Hungary were sam¬ 
pled. Whole blood was obtained by foreleg 
vein puncture. Blood samples were stored in 
a special EDTA -Thymol buffer at minus 20 
°C. Captured individuals were marked with a 
dot of nail polish, thereby preventing recap¬ 
ture of the same individuals. Total genomic 


30 


Genetic differentiation and diversity of Lacerta viridis viridis 


Regions 

Pop 

Location 

Sample size 

Germany, Brandenburg (D) 

1 

secured 

6 


2 

secured 

3 


3 

secured 

3 


Cp 

Linum 

5 

Czech Republic, Bohemia (CzB) 

4 

Karlik 

4 

Czech Republic, Moravia (CzM) 

5 

Podny 

2 


6 

Pavlov 

2 


7 

Bzenec 

3 

Austria, Lower Austria (A) 

8 

WeiBenkirchen 

3 


9 

Hundsheim 

2 


10 

Gumpoldskirchen 

11 

Hungary (H) 

11 

Godolld 

4 

Slovakia, Slovak Carst (SK) 

12 

Tuma 

9 


Tab. 1. Sampled regions, populations and individuals (sample size) of the present study. Populations were 
numbered consecutively (numbers occur in Fig. 6) Because of their high protection status (conservation 
guidelines) we can not include the exact geographic information of Brandenburg localities. Brandenburg 
population numbers refer to terminology of Elbing (2000). Cp stands for captive breeding population in 
Germany (Brandenburg, Linum). 


DNA was extracted using the Qiagen Blood 
Kit, following the manufactures protocol. 


Amplification and sequencing 

To investigate the regional differentiation 
within L. v. viridis, we amplified an mtDNA 
fragment encompassing the partial cyto¬ 
chrome b gene, tRNA Thr, tRNA Pro and 
partial control region (CR). This amplifica¬ 
tion was conducted for 47 individuals with 
the primers LvF and LvRl (Fig. 1, Tab. 2) 
under the following conditions: 1 x buffer 
(Sigma), 1.5 mM MgCl 0 , 0.2 mM dNTPs, 15 
pmol of each primer and 1 U Taq polymerase 


(Sigma). The PCR was performed on an Ep- 
pendorf Mastercycler for 35 cycles with 95 
°C for 30 s, 60 °C for 30 s, 72 °C for 3 min. 
For sequencing, four internal primers (Tab. 2) 
were designed. This resulted in a fragment of 
a final length of about 2600 bp. Purified PCR 
products were used for a cycle sequencing 
reaction using Terminator Ready Reaction 
Mix 'Big Dye' Version 3.1 (Applied Biosy¬ 
stems) following the manufacturer’s protocol 
and were analysed with an ABI 3100 DNA 
Sequencer. The connection of partial se¬ 
quences and internal sequences was done by 
using the program DNASIS v.7.0 (Hitachi 
Software). Final sequences were checked by 
visual inspection using the program Bioedit 


LvF LvFl SP 12Lint4 LvRl 


► ► 




cytochrom b (851 bp) 

1 

P 

• 1 control region (1547-1687 bp) 


(2538 - 2680 bp) 

Fig. 1. Composition of the amplified (primers LvF and LvRl) mitochondrial fragment including the partial 
cytochrome b gene, tRNA Threonine (T), 17 bp Spacer (SP), tRNA Proline (P) and a partial fragment 
of the control region (CR). The dotted area at the 5 'end of the CR comprises a length variable repeat region 
consisting of a 35 bp repeat unit, amplified with primers LvFl and 12Lint4. Internal Primers (Tab. 2) are 
not shown. 
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Primer Sequence 

LvF 5 - CTGCATTTACCTCCATATTGGACG - 3' 

LvFl 5'- GCCTATGCAATCCTTCGCTC - 3' 

LvF2 5'- CAGTTATGCTATGAGCAAGGGTA - 3' 
LvRl 5'- GGCTTTGTAGTTTRATCCTGAC - 3' 

12Lint4 5'- TACCCTTGCTCATAGCATAACTG - 3' 

LvR3 5'- GGGCGGAATGTTAAGGTCCGTTG - 3' 


Tab. 2. Primer used in this study to sequence the mtDNA regions including the cytochrome b gene and 
mitochondrial control region (CR). The Primer 12Lint4 was designed by Brehm et al. (2002), LvF is a 
modified version of cBL (Brehm et al. 2002). The remaining primers were newly designed for this study. 


Version 7.1 (Hall 1999). GenBank accession 
numbers for nucleotide sequences reported 
here are AM087227, AM 087228 and 
AM087289 - AM087330. Final sequence 
alignment was done using the program Clus- 
tal X (Thompson et al. 1997). 


Phylogenetic relationships throughout the 
sampled range 

To analyse the phylogenetic relationships 
and to test whether the populations of L. v. 
viridis in Brandenburg could represent a dis¬ 
tinct subspecies as suggested by Hecht 
(1930), we aligned ten individual L. v. viridis 
mtDNA sequences from all over the sample 
range, one individual of Lcicerta viridis meri- 
dionalis (AM087227) from Greece and one 
individual of Lacerta agilis (AM087228) 
from Romania. Lacerta agilis was used as the 
outgroup. For the construction of a Maxi¬ 
mum-Likelihood phylogenetic tree (ML) we 
used PAUP* version 4.0 blO (Swofford 
2002). DNA substitution rate was calculated 
with a hierarchical likelihood ratio test using 
Modeltest 3.6 (Posada & Crandall 1998). 
The selected model was HKY +I+G (Hasegawa 
et al. 1985) with base frequencies of A = 
0.2946, C = 0.2544, G = 0.1153, T = 0.3357; 
proportion of invariable sites of I = 0.6708 
and a gamma distribution shape parameter of 
variable sites of G = 0.4713. A Neighbor- 
Joining tree (NJ) using Tamura-Nei model 
(Tamura & Nei 1993) and a Maximum-Parsi¬ 
mony tree (MP) with heuristic search using 
ten stepwise additions of sequences and TBR 


branch swapping option was also constructed 
using PAUP* 4.0 blO. To test the robustness 
of NJ and MP bifurcations, bootstrap analy¬ 
ses with 2000 replicates were performed. 


Genetic differentiation between regions 

To estimate the relative migration rates and 
genetic differentiation between the regions, 
all aligned individual sequences were as¬ 
signed to sequence groups representing the 
six main regions (Tab.l). Geneflow and ge¬ 
netic distances between regions were calcu¬ 
lated using DnaSP 4.0 (Rozas et al. 2003). 
Significance of population differentiation 
was tested using Chi-square (Nei 1987, Hud¬ 
son et al. 1992) and permutation test (Hudson 
et al. 1992) with 999 iterations. To obtain a 
pattern of the differentiation between the 
mtDNA haplotypes, the matrix of the pair¬ 
wise genetic distances (F ) between regions 
was analysed using a Neighbor-Joining algo¬ 
rithm (Saitou & Nei 1987) utilising MEGA3 
(Kumar et al. 2004). Within this analysis, 
gaps caused by the differing repeat numbers 
near the CR were excluded. 


Intra-populational variation 

To estimate the genetic diversity within po¬ 
pulations we focussed on the most variable 
part of the analysed mtDNA fragment: the 
repeat region of the CR (Fig. 1). The haplo- 
type diversity of a given population corre¬ 
sponded therefore to the number of different 
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repeats. To increase the data set of this sen¬ 
sitive parameter, we sequenced ten addi¬ 
tional individuals for the CR repeat region 
with the internal primers LvFl and LvR2 
(Fig. 1). The rest of the protocol was identical 
to the analysis of the population differentia¬ 
tion. The analysis comprised twelve different 
populations and one captive population 
within the six sampled regions, which were 
described above (Tab. 1). 


Results 

The general organisation of the amplified 
fragment displays the typical vertebrate pat¬ 
tern (Fig. 1, Boore 1999) consisting of 851 bp 
of cyt b, tRNA Thr (67 bp), tRNA Pro (69 bp) 
and partial control region (CR, 1528-1670 
bp). Both tRNA's are separated from each 
other by the insertion (Fig. 1, spacer) of 
additional seventeen nucleotides. These nu¬ 
cleotides have not been observed in other 
squamate mtDNA sequences available in the 
GenBank e.g. Lacerta dugesii (Brehm et al. 
2003), Podarcis siculus (Podnar et al. unpub¬ 
lished), Pariocela eregia lividus (Kumazawa 
& Nishida 1999) and Cordylus warreni (Ku¬ 
mazawa 2004). Therefore, an assignment of 
these nucleotides to either tRNA Thr or 
tRNA Pro was impossible. We observed an 
overall decreased G content in the light 
strand nucleotide frequencies (A 29.2- 
29.8 %; G 11.5-12.1 %; T 33.1-33.8 % and C 
25.1-25.7 %), which is consistent with the 
known reptilian pattern so far (Macey et al. 
1997). The organisation of the CR is very 
similar to the control region of L. dugesii 
(Brehm et al. 2003). However, the mitochon¬ 
drial genome of L. v. viridis incorporates a 
distinctive and remarkable character: a re¬ 
gion consisting of variable number of repeat 
units next to the tRNA Pro at the 5 'end of CR 
(Fig. 1). The 35 bp motive of this repeat unit 
is highly conserved within L. v. viridis , but 
differs in repeat number (6-10). This size 
polymorphism in form of a variable number 
of tandem repeats (VNTR’s, Lunt et al. 1998) 
results in a varying overall fragment length 


from 2538 bp to 2680 bp between individu¬ 
als. Different types of mtDNA molecules 
within one individual (heteroplasmy) shown 
in other studies and reviewed in Rand (2001) 
were not observed in L. v. viridis so far. The 
mentioned VNTR’s are often observed in 
mitochondrial CR (Sumida et al. 2000, Mundy 
& Helbig 2004) and are supposed to yield 
information about population differentia¬ 
tion and biogeography (Lunt et al. 1998). 


Phylogenetic relationships throughout the 
sampled range 

All phylogenetic trees (NJ, MP and ML) 
inferred from aligned complete mtDNA se¬ 
quences were summarized in Fig. 2. Lacerta 
v. viridis form a monophyletic group distinct 
from L. v. meridionalis, which is supported 
by bootstrap values of 79 % for NJ, 72 % for 
MP and 62 % for ML (Fig. 2). This implicates 
that all sampled regions are inhabited by the 
same subspecies L. v. viridis and that the 
northernmost populations in Germany 
(Brandenburg) do not represent a distinct 
subspecies. 


Genetic variation and differentiation of 
regions 

Altogether, the genetic variation of the ana¬ 
lysed fragment is low. Of 2680 aligned sites 
only 79 were variable and only 44 sites of 
these yielded a phylogenetic signal 
(MEGA3). Generally, the observed differen¬ 
tiation pattern in Fig. 3 fits well to the geo¬ 
graphic distribution of the Green Lizards in 
the studied regions. Pairwise genetic dis¬ 
tances F (permutation test 999 replicates 
p = 0,000; X 2 = 0,0003, df = 150) between the 
regions correlate significantly with their 
geographic distance and Fig. 4 illustrates 
this positive relationship of genetic and geo¬ 
graphic distance between regions. The re¬ 
gions of Brandenburg (eastern Germany) and 
Bohemia (northern Czech Republic) dis¬ 
played the lowest genetic distance (Tab. 3) 
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0.01 




Lacerta viridis 


Fig. 2. Neighbor Joining tree (PAUP* 4.0 blO) obtained from 10 L. v. viridis mtDNA (cyt b, CR) 
sequences within the sample range (D- Germany (Brandenburg), CZ- Czech Republic, SK- Slovakia, H- 
Hungary, A- Austria), one L. v. meridionalis and the outgroup Lacerta agilis. Distances shown are 
Neighbor-Joining distances. Numbers (NJ/ MP/ML) at the nodes show bootstrap values for NJ analysis 
(first number) bootstrap values for Maximum Parsimony (second number) and Maximum Likelihood 
probabilities (> 50 %) resulting from 2000 replications. 


although they are not the geographically 
closest regions within the sampled range. 
Interestingly, the geographically closer re¬ 
gions of Austria and Moravia (southern 
Czech Republic) as well as Slovakia and 
Hungary are genetically more distant to each 
other (Fig. 3). 


Genetic diversity within populations 

Despite an observed low overall nucleotide 
variation within the investigated mtDNA 
fragment, there appear remarkable variations 
within the repeat region of the mitochondrial 


control region (CR, Bohme et al., unpub¬ 
lished). However, the variation was mostly 
limited to insertions or deletions of repeat 
units. This high mutation rate can be traced 
back to a special D-loop structure and 
slipped-strand mispairing during mtDNA re¬ 
plication (Hoelzel et al. 1994, Sbisa et al. 
1997, Savolainen et al. 2000). Changing re¬ 
peat number was found across the whole 
sampled species range but also within popu¬ 
lations. An analysis of the variable CR repeat 
number and the observed haplotype frequen¬ 
cy showed five different repeat haplotypes 
within the investigated populations. Further¬ 
more, we detected that smaller repeat num- 
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D (1-3) 



Fig. 3. Neighbor Joining (NJ) distance tree based 
on pairwise genetic distances (F ST ) of the mtDNA 
haplotypes between six regions within the range of 
Lacerta v. viridis. Letters refer to regions and the 
numbers in parenthesis refer to numbered popula¬ 
tions (Tab. 1, Fig. 6). 



Fig. 4. Correlation of pairwise genetic distance 
(F st ) and geographic distance between different 
regions within the range of Lacerta v. viridis. 


bers (r6-r7) show a higher overall haplotype 
frequency (Fig. 5) than larger repeat numbers 
(r8-rl0). Within populations, considerable 
differences occur in the repeat haplotype 
distribution (Fig. 6). More centrally located 
populations, like population 9 (East Austria) 
and 12 (Slovakia), contain private repeat 
haplotypes, haplotypes which only occur in 
one population, e.g. r9 and rlO. Interestingly, 
the most northern populations 1-3 (Ger¬ 
many) and 4 (Czech Republic) contain no 
private repeat haplotypes. 


Region 

D 

CzB 

CzM 

H 

SK A 

D 


0.106 

0.251 

0.444 

0.376 0.416 

CzB 

226 


0.327 

0.496 

0.433 0.461 

CzM 

395 

250 


0.288 

0.126 0.200 

H 

608 

483 

247 


0.208 0.226 

SK 

590 

554 

334 

183 

0.365 

A 

458 

258 

102 

248 

374 


Tab. 3. Genetic differentiation between different 
regions of the northern Lacerta v. viridis subspe¬ 
cies’ range. Values in the upper part display pair¬ 
wise genetic distances (F ST ) between regions, 
whereas the values in the lower part display appro¬ 
ximate geographic distance in km between regions. 
Regions: D - Germany (Brandenburg), CzB - 
Czech Republic (Bohemia), CzM - Czech Republic 
(Moravia), FI - Hungary, SK - Slovakia, A - 
Austria. 


Situation in Brandenburg 

To obtain more information about the endan¬ 
gered and isolated populations in Branden¬ 
burg, twelve individuals from three wild po¬ 
pulations (Popl, Pop2, Pop3) and five indi¬ 
viduals of the captive population were 
analysed. Within the 2680 aligned base pairs 
we found 22 variable sites between the indi¬ 
viduals of the wild populations; of these, 
sixteen yielded phylogenetic information. 
Within the captive population, we found 
only four variable sites with no phylogenetic 
information. In comparison to the wild popu¬ 
lations in Brandenburg (Fig. 7) the captive 
population displayed the strongest relation¬ 
ship to the population Pop3 (F SJ 0,067), 
whereas the genetic distance to Popl and 
Pop2 is higher (F SJ 0,241 / 0,584). Overall 
haplotype frequency of the repeat number in 
Brandenburg populations was low (Fig. 6); 
only two different repeat haplotypes were 
observed (r6 and r8). Like in most of the 
sampled Brandenburg populations, only one 
repeat haplotype (r6) occurred in the captive 
population. 


Discussion 


A genetic analysis of all mtDNA haplotypes 
revealed no strong nucleotide differences 


35 








Manja U. Bohme et al. 



number of repeats 

Fig. 5. Overall haplotype frequency of repeat num¬ 
bers within the mitochondrial control region (CR) 
of Lacerta v. viridis. 

between the northern edge or isolated popu¬ 
lations and the central populations in the 
observed subspecies’ range (Fig. 2) of Lacer¬ 
ta v. viridis. Our results of a positive correla¬ 
tion between genetic and geographic dis¬ 
tance (Fig. 4) suggest a stepwise dispersal of 
individuals over short distances. This fits to 
the hypothetical postglacial migration be¬ 
haviour of L. v. viridis (Noll 1878, Peters 
1970), with a strong relationship to microcli¬ 


matic favourable regions, following rivers 
and the southern slopes of hills. Neverthe¬ 
less, the genetic differentiation between 
Brandenburg (eastern Germany) and Bohe¬ 
mia (Czech Republic) needs to be discussed 
in detail (Fig. 3). There is a remarkably small 
genetic distance between Bohemia (CzB) 
and Brandenburg (DD) in contrast to the 
geographically closer related region pairs 
(Tab. 3) such as Austria (A) and Moravia 
(CzM) or Slovakia (SK) and Flungary (H). 
This result seems to yield information for the 
reconstruction of postglacial migration 
routes. The common hypothesis (Peters 
1970) assumes a Holocene immigration of L. 
v. viridis from the Czech Republic (Bohemia 
or Moravia) to eastern Germany. Our results 
support this hypothesis but also favour a 
route from Bohemia, along the valleys of 
rivers Vltava (Moldau) and Labe (Elbe). 

Generally, in contrast to other vertebrate 
studies which do not focus at the subspecies 
level (Gubitz et al. 2000, Paulo et al. 2002, 



Fig. 6. Haplotype pattern of the variable number of repeats (VNTR’s) within populations of Lacerta v. 
viridis in different regions of its northern range. Circles display haplotypic frequencies within the sampled 
populations. Population numbers refer to Table 1. Populations 1-3 Germany (Brandenburg), 4 Czech 
Republic (Bohemia), 5-7 Czech Republic (Moravia), 8-10 Austria, 11 Hungary, 12 Slovakia. 
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Captive 



Fig. 7. NJ tree based on pairwise genetic distances 
(F st ) of mtDNA haplotypes between Czech (Bohe¬ 
mian) and German wild populations of Lacerta v. 
viridis and the captive population in Brandenburg. 

Brehm et al. 2003, Hirota et al. 2004), the 
mitochondrial DNA (mtDNA) fragment in¬ 
cluding the cytochrome b gene and the con¬ 
trol region (CR) displayed a low genetic 
variability within L. v. viridis. It yields not 
enough phylogenetic information to obtain 
a phylogeny based on individual haplo¬ 
types. The genetic variation within popula¬ 
tions, analysed using the repeat haplotype 
frequency of the CR (Fig. 6), does not show 
a higher genetic variation within the central 
populations in contrast to the northern, more 
isolated populations. There is an overall 
small variation of only five different CR 
repeat haplotypes (Fig. 5). Shorter repeat 
numbers (r6/r7) seem to be the main repeat 
haplotypes whereas longer repeat numbers 
are more uncommon. Furthermore, we de¬ 
tected private repeat haplotypes only in the 
more central populations for example in 
Pop9 and Popl2. Because of the chosen 
amplification conditions we can exclude a 
skewed distribution to smaller size classes 
due to a laboratory artefact as mentioned in 
Lunt et al. (1998). Despite the well known 
high variability of the CR, this marker sys¬ 
tem could not resolve the genetic variation 
of the populations satisfactorily. Further in¬ 


formation of nuclear microsatellite data 
(Bohme et al. 2005) is needed to clarify this 
genetic situation (work in progress). How¬ 
ever it is important to note that due to the 
comparison of the endangered populations 
in Brandenburg with the central populations 
in Slovakia and Hungary, the low variability 
of the CR was revealed. A single analysis 
which would have focussed only at the en¬ 
dangered populations would have led to the 
false conclusion that the genetic diversity of 
these populations is strongly reduced where¬ 
as in reality the molecular marker itself dis¬ 
plays a low variation. Therefore this is an 
important example that genetic markers need 
to be „calibrated“, before they can be used to 
estimate the genetic diversity of endangered 
populations. 

Captive and Brandenburg populations 

We observed an overall low genetic mtDNA 
variation in the captive population and the 
Brandenburg populations and no pro¬ 
nounced differences with other northern po¬ 
pulations within the subspecies range of L. v. 
viridis. Therefore, genetic data (Fig. 2) con¬ 
firm the morphological findings of Mertens 
& Schnurre (1946, 1949), that the Branden¬ 
burg populations represent the subspecies L. 
v. viridis and do not represent a distinct 
subspecies as suggested by Hecht (1930). 
Also the CR repeat data (Fig. 6) support this 
conclusion. 

Within the Brandenburg populations our 
data (Fig. 7) show a strong genetic link be¬ 
tween the captive population and Pop3. This 
is expected, as the founder individuals of the 
captured population originated from this po¬ 
pulation (Kirmse 1990, 1994). Consequent¬ 
ly, we are confident that the genetic methods 
employed by us are appropriate to analyse 
the genetic relationships between popula¬ 
tions in L. v. viridis. There is also a surprising 
genetic differentiation between the three 
wild populations (Popl/2/3), suggesting that 
within a small geographic distance of only 2- 
14 km, genetic differentiation of populations 
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is possible. Maybe the small individual num¬ 
bers, the low population density (Nettmann & 
Rykena 1984, Elbing 2001b) and the strong 
microrelief (wall) in some localities will en¬ 
force this effect. 
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